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Abstract 

Dependence on the temperature of the integration constants of the relativistic 
quantum Hamilton-Jacobi leads to Lorentz symmetry breaking, a consequence of the 
nonlocal nature of the quantum potential. There is an intriguing analogy with the 
spontaneous symmetry breaking of the Lorentz group in the context of thermal QFT 
and KMS states, a basic issue in noncommutative geometry. Interestingly, also in 
thermal QFT the breaking of symmetries, including Lorentz symmetry, occur even in 

cases with no interactions. Here we consider the correction term yjrat? /T:KTe~ 'kt . 
Using the estimate value of 20K/km for the Earth's crust thermal gradient, one gets 
rriu^ « 0.46 eV that reconciles OPERA, MINOS and SN1987A data and explains 
why superluminality is not observed in other experiments. We also show that the 
observed deviation of the speed of light (c^ — c)/c 1.7 x 10~^^ would correspond 
to a photon mass ~ 2.3 x 10~^^ eV. This would provide a natural explanation 
of the apparently paradoxical consequences of violating Lorentz symmetry. Finally, 
we derive the Lagrangian leading to the above expressions for the velocity and 
show that the change of the relativistic equations is simply described by means of 
a thermal coordinate. 
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A possible explanation of OPERA's results [T] has been recently proposed in [2]. There 
it has been observed the curious coincidence between OPERA's results and the kinematics 
defined by the relativistic quantum Hamilton- Jacobi equation, derived by first principles 
in [21 m E]. Such a formulation is strictly related to classical-quantum duality [6]. 

Consider the stationary Klein-Gordon equation {—h'^c^A + m^c"^ — E'^)tlj = 0. In 
one-dimensional space the associated quantum Hamilton- Jacobi equation is 

{d.Sof + m'c' - ^ + ^{So, q} = 0, (1) 

where {f,q} = — K/'")^ denotes the Schwarzian derivative of /. The term Q = 

-^{So,q}, is the quantum potential and p = dgSo the conjugate momentum. It is crucial 
that Q, which is always non-trivial [Sj HJ [5], is quite different with respect to the one by 
de Broglie and Bohm. In ^ it has been shown that, upon averaging the mean speed on 
the period [g, g + f], k = \/E'^ — rn^d^/hc, one gets v = c^^^^^^ / Li with Li < 1 an 
integration constant of ([T]). Dependence of the initial conditions of ([T]) on the temperature 
T leads to contributions to the dispersion relation depending on T . We will show that 



V = c (1 + ^/mc^/nKTe kt ) ^ 

E 

reconciles OPERA, MINOS and SN1987A, giving m^^ ^ 0.46 eV and v = c for massless 
particles, if any, and explains why the Cohen-Glashow argument [8] does not apply to the 
above expression for v. Furthermore, since the next lightest known massive particle, the 
electron, is approximately 10^ heavier than m^, by e~KT = (^e~ "f«^ , one sees that 
to have measurable effects on the electron one would need a temperature T lO^K, so 
explaining why corrections with respect to experimental data on the Earth have not been 
observed. Also note that the maximal speed is VMax = c + c/ \/2Tce, which holds only for 
massive particles and is reached at T = 2mc^/K. Another remarkable fact is that the 
observed deviation of the speed of light (c^ — c)/c ~ 1.7 x 10~^^ [7J would correspond 
to a photon mass ~ 2.3 x 10~^^eV. Even an extremely thin value of the mass 
photon would provide a natural explanation of the apparently paradoxical consequences 
of violating Lorentz symmetry. 

Since the quantum potential plays the central role in our construction, it is worth 
stressing that the derivation and formulation of the quantum version of the Hamilton- 
Jacobi equations of |3l is quite different from the one by de Broglie-Bohm. In particular, 
the quantum potential is basically different. Let us consider the case of non- relativistic 
quantum mechanics, de Broglie-Bohm made the identification t/j = Ren^° with ip the 
wave- function, not just a general complex solution of the Schrodinger equation, as it should 
be on mathematical grounds [3l H] . This leads to the following trouble. Consider the case 
in which the wave-function is the eigenfunction of the Hamiltonian corresponding to a 
one-dimensional bound state, e.g. the harmonic oscillator. In this case ijj is proportional 
to a real function. It follows that Sq is a constant, so that the conjugate momentum 
is trivial. This would imply that at the quantum level the particle is at rest and starts 
moving in the classical limit. This is a well-known paradox observed by Einstein. In the 
derivation considered in (SJH] such a paradox is naturally resolved by construction: {Sq, q} 
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is well-defined if dqSo is never vanishing. In particular, the general solution of the non- 
relativistic quantum stationary Hamilton- Jacobi equation, which is formally equivalent to 
([I]), with (m^c"^ — E'^)/2mc^ replaced hy V — E, is 

2i c sax ,„W + ii 
w — It 

with w = ftp G M, where now and tp^ are two real linearly independent solutions of 
the stationary Schrodinger equation. 6 = {a,i}, a E M. and i = ii + U2 are integration 
constants. Note that ^1 7^ even when E = equivalent to having So 7^ cnst, which is 
a necessary condition to define {So,q\. This implies a non-trivial So, even for a particle 
classically at rest. One has = S'q'^^"^ (^AeT^^° + i?e~s'5"), and G M implies \A\ = \B\, 
so that So is non-trivial. The ^ — )■ limit is subtle and leads to the appearance of funda- 
mental constants [31 U |5]. A basic observation in [2J is that averaging the period of the 
oscillating terms, which also appear in the non-relativistic case, leads to /i-independent 
solutions that, besides including the standard one, describe other solutions depending on 
the values of the integration constants (see also later). We note that the above decompo- 
sition of ip, now called bipolar decomposition, is successfully used in studying molecular 
trajectories [9l [10] (see also [TT]). 

Another distinguished feature of the formulation introduced in (Sj H] is that energy 
quantization follows without any axiomatic interpretation of the wave-function. As a 
result, the quantum potential, intrinsically different from the one by de Broglie-Bohm, is 
always non-trivial. Like mc^, it plays the role of intrinsic energy and it is at the basis of 
the quantum behavior. For example, besides energy quantization, it makes transparent its 
role in the tunnel's effect, where guarantees that the conjugate momentum always takes 
real values. 

The general solution of ([T]) is ([2]), where now ip and ip^ are two real linearly independent 
solutions of the Klein-Gordon equation. Following Floyd [13], time parametrization is 
defined by Jacobi's theorem t = Since p = dgSo, it gives the group velocity v = 

dE/dp. Set L = k£, Li = 3fJL, L2 = '^L. The mean speed is 



q \/ E'^ — m'^c^ I sin(A;g) — iLcos{kq)\'^ 
t ~ ^ E Ti 



that for L = 1 reduces to the classical relativistic relation v = ^ . In the following 

we will consider the case \L\ = 1, so that 

VE^ - m'^c^ 1 + L2 smi2kq) 
" = ^ E • 

The integration constants Li and L2 may depend on particle's quantum numbers, energy 
and fundamental constants as well. As we said, here we will consider an intriguing depen- 
dence of L on the temperature. Such a possibility is related to a new way of considering 
the /i — )■ limit which has been introduced in [2]. Let us first note that because of the 
term in sin(2fcg), that typically is very strongly oscillating, the ^ — )■ limit is not 
well-defined. One possibility, considered in [21 S] , is that the integration constants may 
depend on E, h and other fundamental constants. In this way the term is cured by 
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a suitable dependence on ^ of L, a procedure that leads to consider the Planck length 
(note that C. has the dimension of a length). The new way of getting the classical limit 
considered in [2j is based on the averaging of the oscillations. This is reminiscent of the 
Dirac's averaging of the oscillating part of the free electron's speed ^cha^e~'^'^^^^^H~^ (see 
Dirac's treatment of the free electron in his book). In [2] it has been observed that in the 
case of the OPERA's experiment, the term sin(2A;g) is very rapidly oscillating, with the 
neutrino Compton wavelength being approximately 10^^ times A;~^. This is a stringent 
physical motivation to average v on the interval [g, g + ^]. By ^ (here we denote (v) of 

m by v) 

k Z-^+f , ^E^ - m^c^ 1 
v = - v{q)dq=c — . (4) 

In [2] it has been noted that this may indicate that neutrinos live on a discrete space-time. 

The basic outcome is that now the classical solution depends on Li, whose value, by 
construction, can be also less than one. This leads to deviations of the standard classical 
dispersion relation. Now note that a possible dependence of L on the temperature should 
be not a surprise. Actually, as observed by Dirac in computing the speed of the free 
electron, the uncertainty principle implies that it does not make sense considering scales 
which are much shorter than the Compton wavelength. On other hand, averaging on the 
period may lead to the breaking of the Lorentz group. This is just reflected by possible 
values of Li which are less than the unity. Furthermore, it is natural to expect that 
such an averaging leads to a dependence on the degrees of freedom of the averaged space 
domain, which in general is not the empty space. In doing this we should remind the 
analogy with spontaneous symmetry breaking of Lorentz group in thermal QFT, we will 
shortly discuss later. Also note that the quantum potential is strictly related to the Fisher 
information and Shannon entropy \l2i (see also [llj). It is then natural to expect that also 
here the temperature plays the central role. For such reasons, further discussed below, 
we will consider the case in which L depends on T. 

Remarkably, dependence of L on the temperature is allowed just as a consequence of a 
basic physical property such as linearity of quantum mechanics. In the case of the Klein- 
Gordon equation, as in the case of the Schrodinger equation, such a property is equivalent 
to the invariance of the Schwarzian derivative under Mobius transformations. In this 
respect, note that a basic identity at the basis of |3llH[5] is that {dgSoY can be expressed 

as the difference of two Schwarzian derivatives = ^ (^{e^'^°,g} — {iSo,g}j, that 

forces us introducing the dimensional constant /3, that is h. We also note the appearance 
of the imaginary factor. It follows that Eq. ([T]) is equivalent to the Schwarzian equation 
{e~^^°,q} = jrT^{E'^ — m^c^), which is invariant under a Mobius transformation of e~^^° 
so that we can also assume that a and i depend on the temperature so that Eq.([2]) has 
the form 

f5o{5(r)} _ iajT] ^ + ^^(^) /P-N 

w- ii{T) ■ ^ ^ 

This provides an intriguing possibility in the case of non-relativistic quantum mechanics. 
Actually, possible dependence on the temperature provides a sort of dynamics non de- 
tected by the wave-function and may open a new view on its collapse. As the temperature 
changes, e~^^°^^'^'^^^ moves on the boundary of the Poincare disk (discrete jumps by Mobius 
transformations are also allowed) with {e~^^°^q} remaining invariant. This happens also 
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when {e^^°, q} is considered in the case of the relativistic quantum Hamilton- Jacobi equa- 
tion (II]). In particular, {e"^'^'', q} = j^{E'^ — rn^c^) is invariant under variations of both 
the coordinate and the temperature. 

Eq. ^ may be related to t — T duality and non-commutative geometry (see [2] and 
references therein). In this respect, it is worth recalling de Broglie view [15]. In particular, 
in the framework of relativistic thermodynamics, he considered tm = h/m(? as a particle 
internal time and identified it with an internal temperature KTm = mc?. de Broglie 
also suggested that, in the case of non- relativistic quantum mechanics, the collapse of the 
wave-function is related to a sort of particle's Brownian interaction with the environment. 

Time-temperature duality is a basic feature of QFT, related to the analogy with clas- 
sical statistical mechanics where the inverse temperature plays the role of imaginary time. 
Time and temperature also mix in performing analytic continuation along complex paths 
in the path-integral. Another well-known analogy concerns the transition amplitude for 
a particle for the time it that coincides with the classical partition function for a string 
of length t at /3 = [16j. It is widely believed that such a dualities are deeply related 
to the properties of space-time and should emerge in the string context. 

In [13] it has been considered a string model where time-temperature and classical- 
quantum dualities, emerge naturally. It turns out that the solitonic sector of compactified 
strings have a dual description as quantum statistical partition function on higher di- 
mensional spaces, built in terms of the Jacobian torus of the string worldsheet and of 
the compactified space. More precisely, in [H] it has been shown that in the case of 
compactification on a circle 

e-^^— = tr e-^^ , (6) 

where /3 = 2R^/a\ with R the compactification radius, and the Hamiltonian H is 
Ajj^/27r, with Aj^ the Laplacian on the Jacobian torus Jq of the worldsheet. Eq.dH]) 
is just the direct consequence of the stronger identity H^rn,n = Sm,n'^m,n, that is the set 
{Sm,n\^,n G Z^} coincides with the spectrum of H, a result deeply related to the theory 
of Riemann surfaces [17] . 

Temperature-time duality naturally emerges as a consequence of the complexified 
version of T-duality, a fundamental feature of string theory. By the standard T- 
duality corresponds to the invariance, up to a multiplicative term given by powers of 
(3, of the partition function under inversion of the temperature /3 — > ^. Complexifi- 
cation of /3 has basic motivations which interplay between physics and geometry. Set 
co{A) = tr (^Ae~^^) / tr e~^^ . Using the invariance of the trace under cyclic permutations 
we have 

uj{A{t)B) = u{BA{t + il3)) . (7) 

In such a context the complexification of /3 naturally appears in globally conformal in- 
variant QFT [18]. It is worth noticing that both in [T8j and in the BC system [19] the 
KMS (Kubo-Martin-Schwinger) states [2U] play a central role. In particular, in the limit 
of 0-temperature the KMS states may be used to define the concept of point in noncom- 
mutative space, a basic issue in string theory [2T| . 

A nice feature of the possible temperature dependence of the integration constants of 
is that it does not arise as a consequence of a specific interaction, rather, as shown 
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above, it is just due to the linearity of the Klein-Gordon equation, or, equivalently, of the 
Mobius symmetry of the Schwarzian derivative. 

The above remarks suggest asking whether there exists an L{T) such that the resulting 
expression for the velocity agrees with the experimental data. Let us then go back to the 
KMS states and note that the Fourier transformed form of ([7j) [22] 

FT[u{AB)]{p) = J d'^xe'P^''-y^uj{A{x)B{y)) = e-^P°FT[u{BA)]{~p) , (8) 

shows that FT[u{AB)]{p) and e~'^'^°FT[u{BA)]{—p) cannot be both Lorentz invariant. 
Spontaneous breakdown of Lorentz boost symmetry at T 7^ is nicely seen considering 
the (1,1) component of the 2x2 matrix propagator of a scalar field that at the tree level 

is [22] 

1 p-/3bo| 

FT[uj{vv)]{p) = -— + 2it6{p' - m')- 



p"^ -m? + ie 1 - e-^IPol ' 

which breaks Lorentz u{\MQi,Tip{x)ip{y)]) 7^ 0. 

Even if we are considering a breaking of the Lorentz group from the quantum potential, 
which is a highly nonlocal term, it is clear that the analogy with thermal QFT suggests 
considering in our approach = e~^, as order parameter of the Lorentz symmetry 

breaking. The crucial property of the KMS states is ([7]), or equivalently ([8]), showing that 
time evolution is invariant, upon commutation, under an imaginary shift of the time which 
is inversely proportional to the temperature. In the case of the quantum Hamilton- Jacobi 
equation, time is generated by applying Jacobi's theorem. Then one gets t as a function 
of q and a shift of time corresponds to a shift of q. On the other hand, this corresponds 
to a jump of e^^° on the boundary of the Poincare disk, and may be compensated by 
changing L{T) that, depending on its functional structure, may correspond to a change 
of T. 

In the present formulation the dependence on e~^<T as order parameter of the breaking 
of Lorentz symmetry comes out as an effect due to the quantum potential and it is not 
seen as due to some particle interaction. For such a reason it is not a surprise that a so 
weakly interacting particle, such as the neutrino, may get thermal contributions to the 
dispersion relation. Remarkably, even in this case there is an analogy with the thermal 
breakdown of symmetries in thermal field theory since the spontaneous breakdown of 
symmetries occurs even in cases with no interactions. 

We then look for an expression of Li leading to a speed which is a agreement with the 
known data. We now show that there is a remarkable solution that passes several tests. 



Namely, by (j4]) it follows that choosing l/Li(T) = 1 + W ^^T^^e kt , leads to 



v-c (l + ^_.-w). (9) 

Note that, just as in the case of the temperature in the Jiittner distribution, which is the 
one measured by an observer at the rest with the gas as a whole, even here the temperature 
is the one measured by an observer at the rest with the medium. Let us show that ([9]) 
passes the known possible tests. First, note that knowing the temperature distribution 
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along neutrino's path 7 one may compute its flight time 



, ds , , 
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Estimating the Earth's crust thermal gradient to be 20K/km (see below), and the room 
temperature 293K at CERN and CNGS, we have T{s)/K = 20(r - Vs^ - Is + r^)10-^ + 
293, where r = [P + Ah'^)/8h, with I = 732km the baseline length and h = 11 Akm the 
maximal deep, we get that (f — c)/c ~ 2.37 x 10~^ corresponds to the neutrino mass 

rrij,^ ^ 0.46 eV , 

which is in the range of its estimated value. Eq.([9]) also agrees with MINOS's results [23] . 
Another important test is that in the case of neutrinos from the SN1987A supernova [21] 
one has T ^ during the flight from the supernova to the Earth, so that 



ttKT^ 







that is for the neutrinos coming from the SN1987A supernova Eq.Q just provides f ~ c, as 
observed. Another crucial test for Eq. is that it gives f = c for massless particles, if any. 
Furthermore, it fixes the maximal speed to be VMax = c + c/ A/27re, which is independent 
of the mass but holds only for massive particles. It is reached at Tmux = 2mc^ /K, which 
is twice de Broglie's temperature. There is another interesting test for ([9]). Namely, since 
the next lightest known massive particle, the electron, is approximately 10^ heavier than 
rrii,, by e~KT = (e~"f?T^) , one sees that measurable effects on the electron need a 
temperature T ^ lO^i^. This would explain why corrections with respect to experimental 
data on the Earth have not been observed. Also note that for m^,^ the maximal speed 
c + c/ y/2Tie is reached at TMax = 2 x W^K, while for the electron and proton such a 
velocity is reached at Tmox = lO^^K and T^tax = 2 x 10^'^ respectively, so that the 
effects are measurable only inside supernovae. 

An intriguing possibility that would answer several questions on the formulation of 
special relativity in the case of superluminality is that massless particles do not exist. On 
the other hand, like a particle in S^, a compact Universe may imply that energy, and 
therefore masses, admit only discretized values with a minimal gap. A hint in such a 
direction is the result by Gharibyan who observed that photons with 12.7 GeV energy are 
moving faster than light by 5.1{lA)mm/ s [7]. This seems due to an energy dependence on 
the vacuum refraction index, which implies a temperature dependence. The measurement 
was done at room temperature, i.e. the accelerator vacuum was filled with 293K blackbody 
photons. By the correction (c^ — c)/c ~ 1.7 x 10^^^ would correspond to a photon 



"7 ^. . , . ... 

mass 



^ 2.3 X lO^^^eV , 



-18 , 



which is much lower than the upper bound 10~ eV reported by the Particle Data Group. 

A remarkable feature of is that it also explains why the argument by Cohen and 
Glashow [8J for the bremsstrahlung — >■ z/^ + e+ + e~ (see also [25j) does not apply to 
(Mi. To see this let us first show that the modified relativistic relation is 



E'' = pVfl + ./^e-wV + mV . (11) 
ttKT 
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Replacing {m^c^ — £'^)/2mc^ hy V — (II]) is formally equivalent to the non-relativistic 
quantum stationary Hamilton- Jacobi equation. It follows that, as observed in [21 S], 
Eq.([T]) can be seen as a deformation of the classical relativistic Hamilton- Jacobi equation 
by a "conformal factor". Actually, noting that {5o,g} = — ((9g5o)^{g, s}, s = 5o(g), we 
see that ^ is equivalent to 

= h'UiSo)] + , (12) 

with U{So) = {g, s}/2, the canonical potential introduced in the framework of p-q 
duality [3l H]. Eq. f|T2|) can be expressed in the form E'^ = {dqS^Yc? + m^c^ where 



(13) 



dqj \j \ — /3^(g) with p^{q) = h^{q, s}/2. Integrating 

dq' 



Eq. (fT3|) indicates that in considering the differential structure one should take into account 
quantum effects on space geometry. Eq. (IT^ is equivalent to 

: ^ /' , - , (14) 

J VE^ - m?c^ J - m?c^ ^ ' 

The above investigation leads to a basic result. Namely, comparison with the dispersion 
relation ( fTTj) shows that the classical limit we obtained by averaging on the short period, 
allowing the solution with the crucial Li < 1, is equivalent to the thermal deformation of 
the coordinate 

lim q = q^ = -i= , 

1 + \ f^e KT 



■kKT^ 

so that the relativistic quantum Hamilton- Jacobi equation (P) becomes 

dqi 

which is ( ITTj) . Eq.( IT^ is clearly related to the relativistic classical action whose modifi- 
cation, leading to (ITT]) , is 



2 / QSq \ 2 , 2 4 

E = — — c -|- m c 



L 



It 

with 

1 

It = 



1 - 



We have 

- '""''^ (15) 
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and 




(16) 



so that we get ( ITTl) and by v 



V 




ttKT 



e KT 




which is, by (ITT]) . Eq.dH]). The dispersion relation f llip is quite different from the one 
considered by Cohen and Glashow. In particular, here E'^ — p^c^ does not correspond to 
any effective mass. Another aspect in Cohen and Glashow argument is that they assume 
that the neutrino mass m^^ can be neglected so that one may use the dispersion relation 



However, as we saw, the case m = is a critical one since ([9]) gives v = c and ffTTl) E = pc. 
In other words, even if very tiny, a mass different from zero is crucial to have superluminal 
solutions since the dispersion relations in the case m ^ and m = differ substantially. 

A possible further test of would be an experiment with the neutrino beam passing 
through the Earth's core whose estimated temperature is 5.700K. Interestingly, a tem- 
perature dependence of neutrino speed would allow geophysical investigations by neutrino 
beams. Another check concerns possible speed fluctuations with respect neutrino's de- 
tection times at OPERA and MINOS. The point is that there is a high thermal gradient 
inside the Earth crust. Temperature increases by as 20 — 25K per kilometer in the upper 
part of the crust, but its gradient is smaller in deeper crust. The highest temperature 
crust is reached at the boundary with the underlying mantle where it ranges from about 
4:70K to 670K. Small changes of the non-inertial forces, due to the periodical motion 
Moon-Earth and Sun-Earth, reflect in a change of the beam trajectory in Earth's crust. 
This should provide a variation of the temperature distribution along neutrino's path 7 

and therefore of its flight time ffTUl) . Note that, due to the e~w term, a small change 
of T would provide a measurable effect on v. Also note that increasing the energy of the 
beam increases the local temperature so that v may have a weak dependence on the beam 
energy. Therefore, it would be interesting to know the thermal gradient for the first tract 
where neutrinos are produced and where the highest temperature is approximately 600K. 
Let us also mention that, as happened in other experiments, one should also check pos- 
sible effects due to tidal forces on the experimental apparatus. Finally, note that Eq.(|9]) 
may admit higher order terms, possibly energy dependent. In this respect it may be of 
interest considering p6]. 

We conclude observing that superluminal neutrinos have been considered in other 
works before OPERA's results, in particular in [27] . After OPERA the problem of super- 
luminal neutrinos stimulated a lot of work (see [28] for a partial list). 




E = vp . 
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